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ABSTRACT 

We present a sample of hard X-ray selected candidate black holes (BHs) in 19 dwarf galaxies. 
BH candidates are identified by cross-matching a parent sample of ~ 44, 000 local dwarf galaxies 
(M* < 3 x 10 9 M 0 , z < 0.055) with the Chandra Source Catalog, and subsequently analyzing the 
original X-ray data products for matched sources. Of the 19 dwarf galaxies in our sample, 8 have 
X-ray detections reported here for the first time. We find a total of 43 point-like hard X-ray sources 
with individual luminosities Z^-iokeV 00 10 37 —10 4 ° erg s -1 . Hard X-ray luminosities in this range can 
be attained by stellar-mass X-ray binaries (XRBs), and by massive BHs accreting at low Eddington 
ratio. We place an upper limit of 53% (10/19) on the fraction of galaxies in our sample hosting 
a detectable hard X-ray source consistent with the optical nucleus, although the galaxy center is 
poorly defined in many of our objects. We also find that 42% (8/19) of the galaxies in our sample 
exhibit statistically significant enhanced hard X-ray emission relative to the expected galaxy-wide 
contribution from low-mass and high-mass XRBs, based on the — M*—SFR relation defined 

by more massive and luminous systems. For the majority of these X-ray enhanced dwarf galaxies, 
the excess emission is consistent with (but not necessarily due to) a nuclear X-ray source. Follow¬ 
up observations are necessary to distinguish between stellar-mass XRBs and active galactic nuclei 
powered by more massive BHs. In any case, our results support the notion that X-ray emitting BHs 
in low-mass dwarf galaxies may have had an appreciable impact on reionization in the early Universe. 
Subject headings: galaxies: active — galaxies: dwarf — galaxies: nuclei — X-rays: binaries — X-rays: 
galaxies 


1. INTRODUCTION 

It is now well-established that massive black holes 
(BHs) reside in the nuclei of essentially ever y giant galaxy 


Kormendy & Ho 

2013), 

in dwarf galaxies 

is unki 


mining the occurrence of massive BHs in present-day low- 
mass galaxies and studying their properties is currently 
our best observatio nal probe of the BH seed population in 
the early Universe ( Volonteri|2Q12 Greene|2012 Natara- 


jan|2014|) . At present, dynamical detections of small BHs 
hi dh 


m distant dwarf galaxies are not feasible and so we must 
search for accreting BHs that shine as active galactic nu¬ 
clei (AGN). Studying the X-ray radiation from massive 
BHs in dwarf galaxies is particularly important for mod¬ 
els of reionizatio n of hydrogen at very high redshifts (e.g., 
Volonteri & Gnedin 2009) 


Over the last decade, systematic searches using optical 
spectroscopy from the Sloan Digital Sky Survey (SDSS) 
have revealed hundreds of low-mass AGN. The first such 
studies focused on finding broad-line AGN with virial 


BH masses Mbh ^ 2 x 10 6 M 0 (Greene & Ho 2004 
2007) and narrow-line AGN in gala xies with absolute a - 
band magnitu des M n > —20 mag (Barth et al. 2008). 
More recently, Reines et al. (2013) pushed to even lower 
BH and galaxy masses by searching the spectra of dwarf 
galaxies with stellar masses M* < 3 x 1O 9 M 0 for both 
broad and narrow-line signatures of active BHs. 


1 Department of Astronomy, University of Michigan, 1085 
South University Ave., Ann Arbor, MI 48109, USA; 
reines@umich.edu 

2 Department of Astrophysical Sciences, Princeton University, 
Princeton, NJ 08544, USA 

3 Hubble Fellow 


The optically-selected samples are likely just the tip 
of the iceberg as they are biased toward relatively high 
Eddington ratios, and toward galaxies with little ongoing 
star formation. High-resolution X-ray and radio observa- 


tions offer powerful alternatives for discovering massive 

BHs in low-mass galaxies (e.g., Gallo et al. 

2008 2010 

Reines et al. 2011; Reines & DelTer 2012; 1 

Miller et al. 

2012; Schramm et al. 2013| Reines et al. 

2014L 


candidates in 

dwarf galaxies by making use of a wealth of archival data 
from the Chandra X-ray Observatory. With its low back¬ 
ground and ability to resolve point sources from diffuse 
X-ray emission, Chandra can detect accretion signatures 
from massive BHs (even nearly quiescent on es) t hat may 
not be accessible at other wavebands (Soria et al. 1(2006 

Cl.L .o llonnol ^ llonnnl IV UIU M dAH A 


Gho sh et al.||2Q08( |Zhang et al.||2009l |Gallo et aT"| 2010 
Pellegrini||2010p. 


This work bears technical similarities to searches for 
ultra hmhnous X-ray sources (ULXs) in dwarf galaxies 

(e.g., Swartz et al.|2008 l Kaaret et al.|2011 _ Prestwich et 

al.||2013| |Plotkin et al.||2014| |brorby et al.|2014|). ULXs 
are off-nuclear X-ray sources in excess of the Eddington 
limit for a 10 M 0 BH (~ 10 39 erg s -1 ), the majority of 
which are now thought to be luminous high mass XRBs. 
Luminosities in this regime, however, are also commonly 
produced by massive BHs accreting at low Eddington 
ratios (|Ho||2009). 

This paper is organized as follows. Our sample selec¬ 
tion is outlined in Section|2l The Chandra data reduction 
and photometry is described in Section [3] In Section [4j 
we present our analysis and results. A summary of our 
conclusions and a brief discussion are given in Section [5] 
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2. SAMPLE SELECTION 

We construct our parent sample of dwarf galaxies using 
stellar masses provided by the NASA-Sloan Atla^l(NSA) 
of local (z < 0.055) galaxies. The NSA provides im¬ 
age mosaics, photometry and various galaxy parameters 
based on a re-analysis of optical and ultraviolet observa¬ 
tions from the SDSS and the Galaxy Evolution Explorer 
(GALEX). Galaxy stellar mas ses in the NSA are derived 
from the kcorrect code of Blanton & Roweis (2007). 
Masses are given in units of M Q h~ z and we adopt h = 
0.73. To select dwarf galaxies, we impose a stellar mass 
upper limit of M* < 3 x 1O 9 M 0 , which is approximately 
equal to the stellar mass of the Large Magellanic Cloud 
(LM C). Our mass threshold is identical to that used by 
Reines et al. (2013) in their optical spectroscopic search 
for active massive BHs, and leaves us with 44,594 objects. 

After selecting dwarf galaxies from the NSA, we cross¬ 
match our parent sample to the list of X-ray sources 
within the Chan dra Source Catalog (CSC, Release 1.1; 
Evans et al. 2010). The CSC includes point and com¬ 
pact (^ 30' ) sources detected in the Advanced CCD 
Imaging Spectrometer (ACIS) and High Resolution Cam¬ 
era (HRC) imaging observations from approximately the 
first eight years of Chandra operations, and a source must 
have a flux at least three times larger than its uncertainty 
to be included in the CSC (see Evans et al. 2010). Using a 
match radius of 5" from the nominal center of the galaxy 
as given by the NSA, we cull an initial 61 objects from 
our parent sample with an X-ray source located close to 
the optical galaxy center. 

We then examine the optical properties of the remain¬ 
ing 61 galaxies in detail to determine if each is indeed 
a dwarf galaxy. We check the absolute magnitudes in 
all SDSS bands given in the NSA and redshifts from the 
SDSS spectra (and/or NED when available), as well as 
visually inspect the SDSS images. This results in the 
identification of many interlopers: nearby massive galax¬ 
ies with erroneous mass estimates (these tend to have an 
obviously incorrect magnitude in an SDSS band), por¬ 
tions of nearby galaxies (including Hll regions and star 
clusters), and faint high-redshift galaxies. Interloping 
massive galaxies hosting AGN, in particular, will prefer¬ 
entially correlate with X-rays. 

For the remaining 31 bona fide dwarf galaxies, we ob¬ 
tain and re-analyze the original X-ray data from the 
Chandra archive as described in below. By re¬ 
analyzing the X-ray data, and not simply using the CSC 
catalog products, we are able to tailor the X-ray anal¬ 
ysis toward our specific needs. Some of the advantages 
(which are described in more detail in Q include the 
following: we consider new Chandra observations that 
were taken after the creation of the latest CSC; we ap¬ 
ply the latest Chandra calibration files to each image; we 
customize source and background photometry apertures 
to each individual X-ray source; we align the Chandra 
astrometry to the SDSS reference frame when possible, 
which is important for determining if an X-ray source 
is consistent with the optical center of each galaxy; and 
we are able to search for off-nuclear X-ray sources within 
each galaxy (the initial match to the CSC only isolates 
X-ray sources 5" from the optical nucleus). 

4 http://www.nsatlas.org 


To minimize possible contamination from diffuse X-ray 
emission (e.g., hot gas from star formation) and point- 
like soft X-ray sources (such as BH and neutron star X- 
ray binaries in their thermal dominant states, super-soft 
X-ray sources and/or cataclysmic variables), we exclude 
sources that are not detected in the hard 2-7 keV band. 
As described below, our final sample of dwarfs with hard 
X-ray selected candidate BHs consists of 19 galaxies and 
is given in Table [l] 

3. CHANDRA DATA REDUCTION AND 
PHOTOMETRY 

We obtained observations of 31 dwarf galaxies from the 
Chandra Data Archive, totaling a cumulative exposure 
tim e of 664 ks, and reduc ed the data with ClA0 version 
4.5 (|Fruscione et al.|2006). As many of the galaxies were 
not targeted themselves, their positions are spread over 
the footprint of the detectors with varying distances from 
the nominal aim-points. We reprocessed the data to cre¬ 
ate new level 2 event files and new bad pixel files, taking 
into account if the observations were taken in FAINT or 
VFAINT mode. We then checked for background flares 
and removed time intervals in which the background rate 
was > 3 a above the mean level. For each observation, 
we improved the Chandra astrometry by cross-matching 
the detected X-ray point sources to the SDSS (which 
has absolute astrometry accurate to < O'/I), and applied 
the resulting bore-sight corrections following |Zhao et ah] 
(2005), as follows. First, we create X-ray images, includ¬ 
ing only the 0.3-7.0 keV range, where Chandra is best 
calibrated. Then, we run a wavelet detection algorithm 
on each activated chip using CIA0 wavdetect with a sen¬ 
sitivity threshold corresponding to one expected spuri¬ 
ous source per field (i.e., for a 1024x1024 pixel image, 
this corresponds to a wavdetect significance threshold of 
^10 -6 ). We cross-match the list of X-ray sources with 
optical sources in the SDSS catalog (with r < 23 mag), 
excluding X-ray sources within three half-light radii of 
the target galaxy. If at least two sources were found in 
common, we generated a new aspect solution file for the 
Chandra observation. This resulted in improved astrom¬ 
etry for ^ 80% of the observations and the corrections 
were in the range ^ O'/Ol — l'/Ol. The final astrometric 
solutions are accurate to between 0'/2 and 0'/5. 

After registering the Chandra images to the SDSS, we 
restrict further analysis to the 2-7 keV energy band to 
search for hard X-ray point sources within three half- 
light radii of our target galaxies (as the contribution from 
background X- ray sources is expected to dominate be¬ 
yond this; see Swartz et al. (2008). We run wavdetect 
on each hard image with wavelet scales of 1.0, 1.4, and 
2.0 pixels, using a 4.5 keV exposure map. As before, we 
adjust the wavdetect sensitivity threshold for each im¬ 
age so that we expect at most one spurious detection per 
field. This yielded a total of 60 hard X-ray sources in 
26 galaxies. For the positional accuracy of each source 
centroid, we estimate 95% e rror circles ( p err ) using the 
empirical formula derived by [Hong et al. (2005| see their 
equation 5); this formula depends on number of hard net 
counts detected by wavdetect, and the location of each 
source on the ACIS detector. The majority of the hard 
X-ray sources (86%) have modest positional uncertainties 
<3". The remaining sources are typically located toward 
an edge of the detector, and can have uncertainties up to 
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TABLE 1 

Dwarf Galaxy Sample 


ID 

(1) 

NSAID 

(2) 

SDSS Name 
(3) 

Other Name 

(4) 

log M* 

(5) 

z 

(6) 

N H 

(7) 

Mg 

(8) 

' - 50 

(9) 

1 

130302 

J014446.40+170634.0 

MRK 0361 

9.230 

0.027 

5.03 

-20.0 

2.53 

2 

131489 

J021403.59+275238.0 

NGC 0855 

9.054 

0.002 

6.41 

-17.3 

16.16 

3 

134049 

J050144.00—041718.9 

IC 0399 

9.343 

0.013 

6.46 

-19.3 

4.91 

4 

135806 

J091721.89+415437.9 

UGC 04904 

9.068 

0.006 

1.03 

-17.7 

12.41 

5 

135954 

J093401.99+551427.9 

I Zw 18, MRK 0116 

6.599 

0.003 

1.99 

-14.9 

3.09 

6 

34462 

J103231.89+542403.5 

UGC 05720 

8.934 

0.005 

0.90 

-19.2 

4.51 

7 

36968 

J103410.14+580349.3 

MRK 1434 

7.354 

0.007 

0.59 

-16.5 

1.86 

8 

90225 

J105120.73+324558.9 

NGC 3413 

8.869 

0.002 

1.99 

-17.6 

12.02 

9 

3264 

J115237.19—022809.9 

UGC 06850 

7.812 

0.003 

2.25 

-16.8 

6.39 

10 

41753 

J121326.02+543631.8 


7.607 

0.008 

1.43 

-15.0 

5.16 

11 

140999 

J121539.20+361937.0 

NGC 4214 

9.071 

0.001 

1.49 

-19.0 

43.08 

12 

30830 

J121923.08+054741.4 

VCC 0344 

9.417 

0.007 

1.58 

-17.3 

2.44 

13 

117784 

J122111.29+173819.1 


8.850 

0.007 

2.68 

-17.7 

6.84 

14 

161603 

J122219.80+300347.9 


8.565 

0.002 

1.73 

-15.6 

15.78 

15 

89394 

J122548.86+333248.7 

NGC 4395 

9.102 

0.001 

1.35 

-18.0 

8.00 

16 

103196 

J124516.87+270730.7 


8.928 

0.004 

0.71 

-18.9 

9.51 

17 

142389 

J124911.59+032318.9 

NGC 4701 

9.241 

0.002 

1.91 

-17.9 

16.51 

18 

31028 

J124957.86+051841.0 

NGC 4713 

9.199 

0.002 

1.95 

-18.1 

23.16 

19 

92811 

J132501.28+362613.7 

NGC 5143 

9.342 

0.019 

0.92 

-18.9 

6.11 


Note. — Column 1: identification number assigned in this paper. Column 2: NASA-Sloan Atlas 
identification number. Column 3: SDSS name. Column 4: Alternative names. Column 5: log galaxy 
stellar mass in units of Mq. Column 6: heliocentric redshift. Column 7: Neutral hydrogen column 
density in units of 10 20 cm -2 . Column 8: absolute g —band magnitude corrected for foreground Galactic 
exctinction. Column 9: Petrosian 50% light radius in units of arcseconds. With the exception of columns 
4 and 7, all values are from the NSA and assume h = 0.73. 


~10" in the most extreme cases. 

We then perform aperture photometry on the hard X- 
ray images. We adopt circular apertures with radii cor¬ 
responding to the 90% encircled energy fraction at 4.5 
keV. The aperture sizes varied from 3"0 - 20 / /6, depend¬ 
ing on the off-axis position of each galaxy. In most cases, 
background counts were calculated from an annulus with 
inner and outer radii of 10" and 30" surrounding the 
source centroid position (masking out overlapping point- 
like X-ray sources when present). For sources farther 
from the aim-point with larger apertures, we used back¬ 
ground annuli with inner and outer radii of 20" and 40". 
The number of counts in each background annulus are 
then normalized to the area of each circular source aper¬ 
ture, to estimate the number of expected background 
counts in each aperture. We use the 4.5 keV exposure 
map to estimate effective exposure times at the location 
of each source on the detector. The net counts were then 
divided by 0.9 to correct for the finite aperture size. 

It is possible for wavdetect to return a small number 
of spurious detections, so we further examine each X-ray 
source. Any spurious detections would be caused primar¬ 
ily by statistical fluctuations (i.e., the optics and detector 
characteristics of Chandra are well-constrained and ac¬ 
counted for during the reduction procedure). Since we 
excluded time intervals subject to unusually high back¬ 
ground levels (and limit ourselves to hard X-ray images 
where the background is typically low; see Table [2]), ran¬ 
dom fluctuations are determined by Poisson statistics. 
Using the source and background apertures customized 
for each X-ray source above, we re-assess which hard 
X-ray sources are detected at a statistically meaning¬ 
ful level. To retain a source, we require a number of 
counts higher than the background level (within each 
source aperture) at the >95% confidence level. For obser¬ 


vations with low sky background (<0.5 expected counts 
within the sour ce apert ure), we set the detection thresh¬ 
old at 3 counts (Gehrels 1986). For sources with large sky 
backgrounds, we adopt the 95% confidenc e level from the 
Bayesian formalism of Kraft et al. (1991), where a typi¬ 
cal threshold is 3-4 countsTor cases with 0.5-1 expected 
background counts within the source aperture. 

Sources below the detection threshold calculated above 
are excluded from further analysis. This resulted in 17 of 
the 60 detected sources being excluded, and the remain¬ 
ing 43 detected X-ray sources are found in a total of 19 
galaxies. The coordinates of each X-ray source are listed 
in Table [2] in order of increasing distance from the nom¬ 
inal optical position of their host galaxy, and we also 
note in Table [2] which sources are also included in the 
CSC. Given our adopted 95% detection confidence limit, 
we expect 5% of the sources in Table 1 2/43) could 

still be false positives due to random statistical fluctua¬ 
tions. Although, 5% is likely an upper limit, given that 
many have counts well-above the 95% confidence limit, 
all were deemed detections by wavdetect, and many X- 
ray sources (30/43) are also included in the CSC. 

Next, we estimate unabsorbed 2-10 keV fluxes using 
version 4.6b of the Portable, Interactive Multi-Mission 
Simulator (PIMM^J) , accounting for Gal actic for eground 
absor ptio n toward each galaxy following [Dickey fc Lock " 


man 


(1990), and assuming a power law spectrum with 


photon index T = 1.8, typica l of accret ing massive BHs at 
moderate Eddington ra tios (lHo|2009l), hard-state XRBs 
as we ll as many ULXs (Swartz et al. 2008; Feng & Soria 
2011). Net count rates/unabsorbed fluxes and hmhnosi- 
ties are given in Table [2] These values may be considered 
lower limits as there may be additional absorption intrin- 


5 See http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html. 
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TABLE 2 

Hard X-ray Sources 


ID 

R.A. 

Dec 

Perr 

Net Counts 

Nbg 

Exp. Time 

C 2 —10keV 

log L-,- 1 cikeV 

d 

CXO ID 

(i) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

1 

26.193495 

17.109842 

3.20 

q i i +3.63 

'*’ L 1 - 2.37 

0.2 

13.70 

4.93 

39.9 

1.42 

J 014446 . 4+170634 

2 

33.516998 

27.877712 

1.06 

q qi +3.62 
°- 0± -2.39 

0.1 

1.54 

46.95 

38.7 

7.88 

J021404.0+275239 

3 

75.433619 

- 4.288618 

0.40 

21 77+ 9 - 37 
' ' -7.07 

0.5 

33.77 

12.80 

39.7 

0.25 

J 050144 . 0—041718 


75.433663 

-4.289359 

0.40 

23.99+4® 

0.5 


14.11 

39.7 

2.47 

J050144.0—041721 


75.431942 

-4.291328 

0.41 

I 8 . 43+45 

0.5 


10.84 

39.6 

11.39 

J050143.6—041728 

4 

139.337673 

41.911788 

0.33 

79.47++° 

0.5 

51.64 

33.25 

39.4 

7.61 

J091721.0+415442 


139.343966 

41.909299 

0.34 

56.13++! 

0.5 


23.49 

39.3 

11.48 

J091722.5+415433 


139.336623 

41.909045 

0.31 

190.58++! 

0.5 


79.75 

39.8 

11.49 

J091720.8+415432 


139.331982 

41.910328 

0.79 

2.80++ 

0.5 


1.17 

38.0 

22.06 



139.346643 

41.906837 

0.41 

15-02+'+ 

0.5 


6.29 

38.7 

22.45 

J091723.2+415424 


139.345338 

41.901709 

0.59 

3 qi+ 3 - 96 
0 - y -2.88 

0.5 


1.64 

38.1 

35.55 


5 

143.508191 

55.241235 

0.57 

5 51+ 4 - 17 
°' OJ --3.42 

0.1 

1.28 

94.33 

39.3 

0.50 

J 093401 . 9+551428 

6 

158.133570 

54.400612 

0.40 

17 54+ 8 - 60 

1 1 -ou _ 6 28 

0.2 

16.38 

21.08 

39.2 

0.81 

J103232.0+542402 


158.131360 

54.403201 

0.51 

6 43 +4 ' 42 

°- q:0 -3.86 

0.2 


7.73 

38.8 

9.79 

J103231.5+542411 

7 

158.542998 

58.062962 

9.73 

34 . 70++ 6 

2.0 

9.47 

89.58 

40.1 

3.02 

J 103410 . 1+580346 


158.539409 

58.063876 

9.85 

23 . 57 +;!! 

2.0 


60.84 

39.9 

5.63 


8 

162.836691 

32.765375 

0.51 

15 49 + 8 - 20 
ro.uy _5 gg 

0.1 

1.73 

193.39 

39.6 

3.43 

J105120.8+324555 

9 

178.155613 

-2.468589 

0.53 

14 39 +7-97 

iu.ou _ 5 62 

0.1 

5.06 

61.57 

39.3 

3.81 

J115237.3—022807 

10 

183.360607 

54.609860 

8.95 

91 51 4 ” 17 ' 37 

6.3 

29.61 

66.81 

40.0 

6.35 

J 121326 . 1+543634 

11 

183.920218 

36.327396 

0.65 

O.OO —4 69 

0.3 

26.03 

7.47 

37.7 

17.50 

J121540.8+361939 


183.908986 

36.328876 

0.41 

45.28+!+ 

0.3 


39.57 

38.4 

17.85 

J121538.1+361944 


183.919917 

36.323634 

0.76 

5.38++ 

0.2 


4.70 

37.5 

19.07 



183.909394 

36.322301 

0.33 

235.24+!+ 

0.3 


205.57 

39.1 

20.50 

J121538.2+361921 


183.921454 

36.315876 

0.62 

6 41+ 4-42 
—3.86 

0.3 


5.60 

37.6 

43.16 



183.900242 

36.312919 

0.65 

|T qO + 4.19 
°- <3Z -3.40 

0.2 


4.65 

37.5 

63.39 



183.936452 

36.312645 

0.52 

10 . 86 +'+ 

0.3 


9.49 

37.8 

81.08 

J121544.7+361846 


183.922478 

36.353828 

0.70 

21.89++ 

0.3 


19.13 

38.1 

101.57 

J121541.4+362114 


183.935980 

36.302550 

0.80 

3.10++ 

0.2 


2.71 

37.3 

106.38 


12 

184.846310 

5.794871 

0.80 

6 . 41 ++ 

1.4 

62.97 

2.00 

38.3 

0.39 

J 121923 . 1+054741 


184.847546 

5.796062 

2.00 

4-20+;!! 

1.4 


1.32 

38.1 

6.50 



184.845120 

5.793374 

1.31 

3.04+;°° 

1.4 


0.95 

38.0 

6.57 

J121922.8+054735 

13 

185.296810 

17.636862 

11.25 

14 30 4-7 ’ 95 

i4,t)U - 5 .60 

1.3 

4.07 

69.98 

39.9 

6.61 

J 122111 . 0+173817 

14 

185.580669 

30.061350 

4.76 

o qq+3.62 

o.o«3 _ 2 3g 

0.1 

4.10 

16.16 

38.6 

3.65 

J 122219 . 3+300340 

15 

186.453622 

33.546937 

0.30 

876 . 17 +!'+ 

0.5 

10.58 

1882.84 

40.1 

0.25 

J 122548 . 8+333248 


186.452117 

33.550481 

0.77 

9 oq + 3.64 

2.80-2 34 

0.5 


6.09 

37.7 

13.77 


16 

191.321862 

27.125506 

0.34 

52 . 00 +!+ 

0.2 

15.82 

68.10 

39.4 

0.40 

J 124517 . 2+270731 


191.320909 

27.126584 

0.40 

18.66+'+ 

0.2 


24.44 

39.0 

4.82 



191.320436 

27.126193 

0.39 

20.89+'!! 

0.2 


27.35 

39.0 

5.24 


17 

192.298020 

3.389506 

0.55 

5-43+'+ 

0.1 

9.00 

1.30 

38.3 

3.05 

J124911.5+032322 


192.291310 

3.379165 

0.73 

q qq + 3.63 
O.ZZ — 2 38 

0.1 


7.71 

38.1 

42.54 


18 

192.481030 

5.319315 

0.98 

4.42+'°! 

0.1 

4.40 

2.23 

38.5 

46.22 



192.473880 

5.311602 

0.55 

7 7^ + 4.60 
‘ • ' Z -4.30 

0.1 


3.90 

38.8 

61.98 

J124953.7+051841 

19 

201.255760 

36.437561 

1.40 

3.30+;!! 

0.1 

4.61 

15.87 

40.1 

2.21 

J132501.3+362615 


Note. — Column 1: galaxy ID (see Table 1). Column 2: right ascension of hard X-ray source in units of degrees. Column 3: declination of 
hard X-ray source in units of degrees. Column 4: 95% positional uncertainty in arcseconds. The values do not include the 0.2-0.5 // uncertainties 
from the astrometric correction. Column 5: net counts in the energy range 2-7 keV, after applying a 90% aperture correction to the net counts 
in each source extraction circle. Error bars represent 90% confidence intervals. Column 6 : number of background counts expected within each 
source extraction circle (after applying a 90% aperture correction). Column 7: effective exposure time in kiloseconds at each source position 
on the detector, determined from 4.5 keV exposure maps. Column 8 : 2-10 keV flux in units of 10 -15 erg s -1 cm -2 corrected for Galactic 
absorption. Column 9: log 2-10 keV luminosity in units of erg s -1 corrected for Galactic absorption. Column 10: distance from optical center 
of the galaxy in units of arcseconds. Column 11 : Name of source in CSC catalog, if present. Sources are listed in i ncrea sing distance from 
the center of the galaxy, and sources with positions consistent with the optical nucleus are highlighted in boldface (see j |4.1| ). 









Fig. 1 . — SDSS images of our dwarf galaxy sample. Red ellipses indicate the positions of hard X-ray sources (from wavdetect). Scale 
bars are 5" in length, except for ID 11 which has a 25" scale bar. 
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sic to the sources. All error bars are quoted at the 90% 
confidence level, unless stated otherwise. 

Finally, we assess the expected number of fore¬ 
ground/background sour ces with in 3_ half-light radii of 
each galaxy. Following Plotkin et al. (2014), we use 
the published X-ray source counts fromthe resolved cos - 
mic X-ray background as given by|Moretti et al. (2003). 
Starting with the hard (2-10 keV) X-ray hux Sdim ' which 
corresponds to the dimmest X-ray source in each field, 
we adopt the expression for the cumulative X-ray flux 
distribution given in Equation 2 by Moretti et al. (2003) 
to calculate N(> S'), i.e. the expected number of X-ray 
sources per square degree having an X-ray flux brighter 
than S. Finally, for each galaxy field, N(> S ) is mul¬ 
tiplied by the area enclosed by 3 half-light radii. With 
the exception of one nearby galaxy with a moderately 
deep exposure (ID 11, NGC 4214; = 4.8), the ex¬ 

pected number of background sources, TVbk, varies be¬ 
tween 0.00025 and 0.74 with a median of 0.016. Across 
the sample, we expect a total of 1.4 background contam¬ 
inants (excluding ID 11). 

4. ANALYSIS AND RESULTS 

We detect a total of 43 point-like hard X-ray sources in 
19 dwarf galaxies (Figure JT1) . Approximately 42% of our 
sample of galaxies (8/19) have X-ray detections reported 
here for the first time - we searched the NED databases 
to identify any X-ray sources that are already known to 
be background/foreground objects or otherwise reported 
in the literature. Notes on individual objects are in the 
appendix. 

4.1. Spatial Distribution of X-ray Sources and 
Candidate Nuclear Sources 

In Figure [2j we show the distribution of the locations 
of the hard X-ray sources within each galaxy normal¬ 
ized by the half-light radius. Roughly half of the X-ray 
sources (23/43) fall within one half-light radius of the 
optical galaxy center (as given by the NS A). Such a dis¬ 
tribution is expe cted for the ( off-nuclear) luminous XRB 
population (e.g., S wartz et al.|||2Qll ). It is possible that 
some of the X-ray detections close to the optical center 
could be AGN and not XRBs. However, since the po¬ 
sitional accuracy degrades for sources falling toward the 
edge of the chip, it is difficult to apply a uniform criteria 
across the entire sample to identify likely nuclear X-ray 
sources. Instead, we determine which X-ray sources have 
95% error circles that overlap with the optical center (we 
add ^0.5" uncertainty from the astrometric correction to 
the Perr estimates), and we find that 11 X-ray sources (in 
10 galaxies) are consistent with the nucleus of the host 
galaxy. We caution that the optical center of a galaxy 
may not be well-defined for irregular galaxies or those 
with star-forming clumps. 

We stress that we consider 11 “nuclear” X-ray sources 
to be a generous upper limit, as many of these sources 
are consistent with the optical center primarily because 
the error circles on their positions are large (mostly 
due to galaxies falling toward an edge of the detector). 
For example, 5 X-ray sources (in galaxy IDs 7,10,13,14) 
have positional uncertainties p e rr > 5" (with X-ray po¬ 
sitions ^3-7" from the optical center), and two of these 
sources fall within the same galaxy (ID 7). These X- 
ray sources could very likely be off-nuclear XRBs. One 
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Fig. 2.— Distribution of the distances of hard X-ray sources from 
the optical center of their host galaxy, normalized by the half-light 
radius. 
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Fig. 3. — Specific star formation rate versus r — i optical color for 
the 19 galaxies in our sample. The color bar indicates r — i color. 

previously unidentified X-ray source (in galaxy ID 1; see 
Appendix) lies 1.4 ± 3.2" from the optical center and 
could be an intriguing nuclear source. The remaining 
sources (IDs 3,5,12,15,16) have much more accurate po¬ 
sitions (p err < 1") and lie <0.5" from the optical cen¬ 
ter. Two of these are newly-identified objects, and three 
are previously identified X-ray sources including the low- 
mass Seyfert galaxy NGC 4395 (see Appendix). 

4.2. Host Galaxy Properties 

We obtained galaxy stellar masses and half-light radii 
from the NS A. The stellar masses span a range of M* ~ 
10 7 —10 9-5 Mq, with a median of ~ 10 9 M©. The galaxies 
in our sample are physically small in addition to being 
low mass, with half-light radii < 2 kpc. The galaxies 
span a range of optical colors and our sample includes 
spheroids, disks and irregular systems (see Figure [l]). 
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TABLE 3 

Expected Galaxy-Wide Luminosity from X-ray Binaries 


ID 

(1) 

FUV 

(2) 

L(FUV) obs 

(3) 

w 22 

(4) 

L(25 ^tm) 
(5) 

log SFR 
(6) 

log L™® keV 

(7) 

1 

16.44 

43.45 

5.01 

43.17 

0.58 

39.80 

2 

15.20 

41.67 

5.39 

40.74 

-1.52 

38.18 

3 

15.42 

43.24 

6.01 

42.15 

0.01 

39.27 

4 

16.54 

42.04 

7.58 

40.76 

-1.23 

38.30 

5 

15.79 

41.64 

7.84 

39.97 

-1.67 

37.54 

6 

14.24 

42.83 

2.53 

42.66 

0.04 

39.26 

7 

16.42 

42.34 

7.45 

41.07 

-0.93 

38.29 

8 

15.10 

41.77 

5.53 

40.75 

-1.44 

38.10 

9 

15.05 

42.22 

4.68 

41.52 

-0.88 

38.34 

10 

19.57 

41.13 

> 9.28 

< 40.40 

-2.22 

37.13 

11 

11.57 

42.51 

2.62 

41.23 

-0.76 

38.59 

12 

21.14 

40.37 

> 8.86 

< 40.42 

-2.98 

38.38 

13 

16.31 

42.32 

> 8.07 

< 40.77 

-1.03 

38.34 

14 

19.73 

40.00 

8.84 

39.50 

-3.00 

37.54 

15 

13.65 

41.75 

5.52 

40.15 

-1.56 

38.20 

16 

13.48 

42.89 

3.96 

41.84 

-0.33 

38.92 

17 

14.73 

42.03 

4.54 

41.25 

-1.10 

38.46 

18 

14.22 

42.14 

5.91 

40.61 

-1.16 

38.41 

19 

17.37 

42.79 

8.12 

41.64 

-0.45 

38.89 


Note. — Column 1: galaxy ID (see Table 1). Column 2: GALEX far-UV 
AB magnitude. Column 3: Observed far-UV luminosity in units of erg s -1 . 
Column 4: WISE 22 |Um Vega magnitude. Column 5: 25 \i m luminosity 
in units of e rg s —1 , using 22 /im flux density (in Jy) as a proxy for 25 /im 
flux density (Jarrett et al.12013]). Column 6: log SFR in units of M© yr -1 . 
Column 7: Total 2-ID keV luminosity expected from low-mass and high-mass 
X-ray binaries. Units are in erg s —1 . 


We estimate star formation rates (SFRs) of the host 
galaxies using a combination of far-UV data from 
GALEX and 22 fi m dat a from the Wide-fie ld Infrared 
Survey Explorer (WISE] |Wright et al.|[2010). GALEX 
FUV magnitudes were obtained from the INS A and in¬ 
frared magnitudes were obtained from the WISE All- 
Sky Source Cat alog (see Table [3| ). F ollowing |Kennicutt 
& Evans | (|2012 ) and Hao et al.| (JzOll ), we estimate dust- 
corrected SFRs as log SFR = log X(FUV) corr — 43.35, 
where L( FUV) cor r = I/(FUV) 0 bs±3.89 L(25 /i m) and the 
luminosities are in units of erg s -1 . While the Hao et al. 


Function (XLF) of XRBs in external galaxies (Fabbiano 


(Grimm et al. 2002 

, 2003 Gilfanov et al. 2004 Mineo et 

al. 2012) and with total stellar mass lor early-type galax- 

ies (iGillanov 2004; 

Humphrey & Buote 

2008! Lehmer et 


bright high-mass XRBs and the latter by longer-lived 
low-mass XRBs. 

Given our heterogeneous sample of galaxies, we esti¬ 
mate the expected collective luminosity from XRB s in 


a given galaxy usi ng stellar masses and SFRs from §4.2 


(poll) relation employs 25 /am luminosities from the In- and the relation in Lehmer _et al. (2010), which accounts 


frared Astronomical Satellite (IRAS), the vast majority 
of our sample is not detected by IRAS at 25 /im. WISE 
is sufficiently more sensitive and all but 3 of our sample 
galaxies have 22 /im detections. Given this, and the fact 
that we expect the flux density ratio between 22 /im and 
25 /im (in Jy)_to be ~ 1 for both late-type and early-type 
galaxies ( Jarrett et al.|2013 ), we use 22 /im flux densities 
as a proxy for 25 /im flux densities. For the 3 galaxies not 
detected at 22 /im in WISE, we derive SFRs using the 
observed (uncorrected) FUV luminosities. The derived 
SFRs are given in Table [3j We plot specific SFRs (SSFR 
= SFR/M*) versus r — % optical color for our sample of 
dwarf galaxies in Figure [3| As expected, bluer (redder) 
galaxies have higher (lower) SSFRs. 

4.3. Expected Contribution from X-ray Binaries 

The luminosities of individual hard X-ray sources in 
our sample are in the range ±2-iokeV ^ 10 37 - 10 4 ° 
erg s -1 . While we are primarily interested in finding 
candidate massive BHs, these luminosities could also be 
attained by stellar mass BHs (or neutron stars) in X- 
ray binaries (XRBs) with hard X-ray spectra (see, e.g., 
McClintock & Remillard 2006). The X-ray Luminosity 


for both high-mass and low-mass XRBs: L ^ keV = 

<aM*±/3SFR, where a = (9.05T0.37) x 10 28 erg s -1 M© 1 
and (3 = (1.62 =b 0.22) x 10 39 erg s -1 (M© yr -1 ) -1 . 
While the Lehmer et al. (2010) relation was derived from 
more massive and more luminous galaxies, the range 
of SSFRs in that work is comparable to our sample 
(SFR/M*. ~ 10 -12 — 10 -8,5 yr -1 ). In Figure [4j we plot 
the total observed hard X-ray luminosity versus the ex¬ 
pected galaxy-wide luminosity from XRBs. 

We caution that there are large uncertainties associ¬ 
ated with th e exp ected total luminosities from XRBs. 


The Lehmer et al. (2010) relation itself has a scatter of 
^0.34 dexT There are also large uncertainties in the stel¬ 
lar masses and SFRs. Uncertainties in the stellar masses 


are expected to be_ 0.3 dex ([Conroy et al.||2009| ) and 
the scatter in the|Hao et al. (201 1|) relation is 0.13 dex. 
Stochastic effects are also a concern for galaxies with low 
stellar masses and/or low SFRs since a small number of 
luminous XRBs can dominate the total X-ray luminos¬ 
ity. The la error bars (±0.4 dex) in Figure [4] account 
for systematic uncertainties in stella r masses and SFRs, 
scatter in the Lehmer et al. (2010) relation, as well as 
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log (expected 2-10 keV luminosity from XRBs [erg s ']) 


Fig. 4.— Left: Total observed 2-10 keV luminosity versus expected total luminosity from low-mass and high-mass X-ray binaries. Points 
are color-coded by r — i color. The solid line shows the one-to-one relation and the dashed line shows ou r med ian X-ray sensitivity limit. 
The lcr error bar accounts for systemat ic uncertainties in stellar masses and SFRs, scatter in the Lehmer et al.|(2010) relation, and errors 
on a and /3 in the [Lehmer et al. (2010} relation. Galaxies with statistically significant excess hard X-ray emission (> 3a, or 1.2 dex) are 
marked with an ad ditio nal circle. Hipfit: Same as left panel, except with X-ray emission removed from point sources consistent with the 
optical center (see j |4.1| ). 


errors on a and /3 in Lehmer et al. (2010). 


4.4. Galaxies with Enhanced X-ray Emission 

Eight galaxies in our sample (42%) exhibit excess hard 
X-ray emission relative to the expected contribution from 
XRBs (Figure |4|. These galaxies have observed hard X- 
ray emission > 3 a (or more than 1.2 dex) higher than 
expected given their stellar masses and SFRs. The X- 
ray sources responsible for the excess emission in five of 
these cases also have positions consistent with the opti¬ 
cal centers of the host galaxies as described in j |4.1| (IDs 
5,7,10,13,15). ID 15 is the well-s tu died dwarf Seyfert 
galaxy NGC 4395 (|Filippenko & Sargent 1989) and ID 5 
is the extremely metal-poor galaxy 1 Zw 18 with a well- 
studied luminous X-ray source (see Appendix). To the 
best of our knowledge, the X-ray sources in galaxies 7, 
10 and 13 are reported here for the first time. How¬ 
ever, we note that these galaxies have somewhat irregu¬ 
lar morphologies without clearly defined nuclei, and the 
X-ray positions are not well-constrained. An additional 
three galaxies (IDs 4,8,19) exhibit enhanced X-ray emis¬ 
sion from non-nuclear sources (see Figure |4|. 

4.5. X-ray Detection Limits 

Given the heterogeneous nature of our galaxy sam¬ 
ple, each galaxy has a different luminosity threshold to 
which an X-ray source could be detected. To estimate 
the minimum detectable luminosity, we assume 3 hard 
X-ray photons (2-7 keV) must be detected within the 
source aperture (corresponding to a 95% confide nce limit 
if there is no contribution from the background; |Gehrels| 
1986). More photons would be required if there is sig¬ 
nificant contribution from the local background, so we 
consider the assumption for 3 photons to provide a ro¬ 
bust lower limit. Then, assuming the effective exposure 


time for each galaxy, the Galactic foreground absorption, 
the distance to each galaxy, and a power law spectrum 
with photon index T = 1.8, we calculate the minimum 
2-10 keV luminosity (L m in) for which we could detect a 
hard X-ray source within each galaxy. The sensitivity 
limits range from 5 x 10 36 to 1 x 10 40 erg s 1 , with a 
median value of 2.3 x 10 38 erg s -1 . The distribution of 
observed X-ray luminosities is largely shaped by these 
detection limits (see Figure |4|. 


5. SUMMARY AND DISCUSSION 

We have constructed a sample of 43 hard X-ray se¬ 
lected BH candidates in 19 nearby dwarf galaxies with 
stellar masses M* < 3 x 10 9 M©. Eight of the galax¬ 
ies in our sample have X-ray detections reported here 
for the first time. After combining the total hard X-ray 
flux from all point sources within each galaxy, nearly all 
of the galaxies have cumulative hard X-ray luminosities 
that are higher than expected from the combination of 
low-mass and high-mass XRBs (given each galaxy’s stel¬ 
lar mass and SFR). This high-level of X-ray emission 
is consistent with expectations, given the X-ray selected 
nature of our galaxy sample, and the sensitivity limit of 
each Chandra observation. For the majority of galaxies, 
we are typically capable of detecting a hard X-ray source 
with Z/ 2 —iokeV ^ 10 38 erg s -1 . Therefore, our study is 
generally only sensitive to the luminous tail of the XRB 
population, and more massive accreting BHs. 

Interestingly, even after considering the above bias, we 
find that a substantial fraction of the galaxies in our sam¬ 
ple (42%) exhibit statistically meaningful enhanced hard 
X-ray emission, the majority of which have X-ray sources 
consistent with the o ptical n ucleus . Tim d warf S eyfert 
1 galaxy NGC 4395 (Filippenko & Sargent 1989) is in¬ 
cluded in this sub-sample, demonstrating that searching 






























X-ray Selected Black Holes in Dwarf Galaxies 


9 


for enhanced (cumulative) hard X-ray emission can in¬ 
deed be a useful diagnostic for identifying AGN in dwarf 
galaxies. However, enhanced hard X-ray emission is nei¬ 
ther a necessary nor sufficient requirement for an accret¬ 
ing massive BH. Enhanced galaxy-wide X-ray emission 
from off-nuclear X-ray sources (thought to be luminous 
XRBs) has been observed in other studies, especially in 
low-metallicity dwarfs where the enhancement is up to 
a factor of 10 larger than for ^solar-metallicity galaxies 

i-lT 


(Prestwich et al. 2013 


Brorby et al. 


Follow-up XYay observations, as well as observa¬ 
tions at other wavelengths can help distinguish between 
stellar-mass XRBs and more massive BHs in our sam¬ 
ple. High-resolution radio observations offer a promising 
way forward since the ratio of compact radio to hard X- 
ray emission is significantly higher for acc reting massive 


BHs compared to stellar mass BHs (e.g., Merloni et al 


2003). The combination of radio and X-ray observations 


have so far revealed c andidate massive BHs in the dwarf 
galax ies Henize 2-10 (Reines et al. 201 H Reines & Deller| 
|2012| ) and Mrk 709 ( Reines et al.| 2014). Eleven galax- 
ies in our sample are detected in either the FIRST^] or 
NVSSQ surveys, however the angular resolution is not 
sufficient to associate radio emission with any particular 
X-ray source. Mid-IR diagnostics (e.g., with WISE) have 
been used to select luminous AGN with high confidence 
(e.g., Stern et ah]|2012 Assef et ah]|2013 ), however using 
mid-IR colors to identify AGIN in dwarf galaxies is not 
straightforward. For example, star-forming dwarf galax¬ 
ies can have m id-IR colors overlapping those of WISE - 
selected AGN (|Izotov et al.||2014 Jarrett et al.||2011|). 
Regardless of the masses of the BHs in our sample, our 
results have important implications for the impact of X- 
ray-producing BHs in low - mass galaxie s on th e epo ch of 


reionization (e.g., Volonteri & Gnedin 2009 Mirabel et 


arf20TT|). 
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APPENDIX 

NOTES ON INDIVIDUAL SOURCES 


1. Mrk 361 (ObsID 6855; PI Komossa). Belongs the field of view of the Ultra Luminous Infrared Galaxy IIIZw035. To 
the best of our knowledge, this is the first paper reporting on the discovery of an X-ray source in the target dwarf galaxy. 


2. NGC 855 (ObsID 7095; PI Swartz). Detection of a ULX in this dwarf galaxy was first reported by Swartz et al. 
(2008). 


(C 

et al. ( 2013| ). See also Tzanavaris et al ||2014| 


3. IC 399 ( ObsID 94 05; PI Gallagher). X-ra y sources in this group dwarf galaxy were first reported by Desjardins 


4. UGC 4 904 (O bsID 6729; PI Kulkarni). ToO observation to follow up on a core-collaps e supernova in the tar - 
get galaxy (Ofek et al.p0 13). Analysis of the point-like X-ray population was first presented by Desjardins et al. (2014). 

5. I Zw 18, Mrk 116 (ObsID 805: PI B omans). A nuclear X-ray source in this metal poo r dwarf galaxy was firs t 
reported by ~Ott, Walter & Brinks ( 2005| ). See also |Kaaret et al. (2011); Kaaret & Feng (2013); Prestwich et al. (2013); 


Brorby et ah 12014): 


6. UGC 572 0 (ObsID 9519; P I Mas-Hes se). X-ray sources a s sociated with in thi s Lym an a emitting galaxy were first 
reported by |Grier et al.| ( 20TT ). See also |Mineo et al. (2012); Otf-Floranes et al. (2012). 


7. Mrk 1434 (Obsid 3347; PI Barger). Belongs to the field of view of the Lockman Hole-Northwest (Yang et al. 


2004). To the best of our knowledge, this is the first paper reporting on X-ray sources associated the target dwarf galaxy. 


The VLA FIRST Survey: Faint Images of the Radio Sky at 


Twenty-Centimeters; |Becker et al. 

(19951 

7 The NRAO VLA Sky Survey, 

Condon et al. 

(1998 
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8. NGC 3413 (ObsID 7102; PI Swartz). The X-ray source in this galaxy is part of the Liu 2011 catalog (it is not, 
however, included in Sw artz et al.||2008~ nor |Swartz, Tennant fe So ria 2009, as both studies focus on ULXs and thus 
discard nuclear X-ray sources). 

9. UGC 06850 (ObsID 7135; P I Swartz) . A bright X-ray s ource in this group dwarf galaxy was first reported by 
Swartz, Tennant & Soria ( 2009| ). See also |Swartz et al.||2011| 


10. SDSS JJ_21326.02+543631.8 (ObsID 7071; PI Kaaret). Belongs to the field of view of NGC 4194 (Kaaret fc 


Alonso-Herrero 2008). To the best of our knowledge, this is the first paper reporting on a X-ray source associated 
with the target dwarf galaxy. 


11. N GC 42 14 (ObsID 5197; PI Zezas). X-ray sources i n thi s nearb y d warf galaxy w e re firs t reported by Ghosh et al. 
(2006) and Leonidaki, Zezas & Boumis ( |2010 ). See also Liu (2011) and Mineo et al. (2012). 


12. VCC 344 (ObsID 956 9; PI Zesaz). Belongs to the field of view of the nearby galaxy NGC 4261 (Worrall et al. 
2010) (see also |Gultek i n et al.| |2009). To the best of our knowledge, this is the first paper reporting on X-ray sources 


in tne target dwarf galaxy. 

13. SDSS J122111.29+173819.1 (ObsID 8085; PI Treu). Belongs to the field of view of VCC 437 (Gallo et al. 2010). 
To the best of our knowledge, this is the first paper reporting on X-ray sources in the target dwarf galaxy. 


14. SDSS J122219.44+300344.1 (ObsID 7853; PI Mathur). Belongs to the field of view of NGC 4308 (Miller et al. 


2012). To the best of our knowledge, this is the first paper reporting on an X-ray source in the target dwarf galaxy. 


15. NGC 4395 (ObsID 5101; PI Moran). This nearby dwarf spiral galaxy harbors the least luminous Seyfert 1 nucleus 
known. O’Neill et al. (2006) first reported on the X-ray population resolved with Chandra . 

16. SDSS J124517.25+270732.1 (Obsid 13928; PI Reines). Reines et al., in prep. 

17. NGC 4701 (ObsID 7148; PI Soria). |Desro ches & Ho (2009) first reported on a nuclear X-ray source in this nearby 


late type dwarf. Here we report on an additional off-nuclear source. See also Liu (2011). 


18. NGC 4713 
this dwarf. Here 


(ObsID 4019; PI Satyapal). Dudik et al.| (2005| first reported on the lack of a nuclear X-ray source in 
e we report on off-nuclear sources. See also Liu (2011). 


19. NGC 5143 (ObbsID 4055; PI Sambruna). Belongs to the field of view of the FRI 1322+36. To the best of our 
knowledge, this is the first paper reporting on X-ray sources in the target dwarf galaxy. 
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